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ABSTRACT 

The fact that a temperature and an entropy may be associated with horizons in semi- 
classical general relativity has led many to suspect that spacetime has microstructure. If 
this is indeed the case then its description via Riemannian geometry must be regarded as 
an effective theory of the aggregate behavior of some more fundamental degrees of freedom 
that remain unknown, in many ways similar to the treatment of fluid dynamics via the 
Navier-Stokes equations. This led us to ask how a geometric structure may naturally arise 
in thermodynamics or statistical mechanics and what evolution may mean in this context. 
In this article we argue that it is possible to view thermodynamic processes as the evolution 
of a dynamical system, described by a quadratic Lagrangian and a metric on the thermo- 
dynamic configuration space. The Lagrangian is an invariant distance between equilibrium 
thermodynamic states as defined by the metric, which is straightforwardly obtained from a 
complete set of equations of state. 
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I. INTRODUCTION 



Research over the past few decades has shown that when the principles of quantum 
mechanics are combined with a geometric theory of gravity in the presence of horizons, 
spacetime exhibits thermodynamic behavior |jl|, Without reference to the gravitational 
field equations, both a temperature and an entropy can be associated with any horizon in 
the following sense: whenever a horizon is present there is also a law, similar to the first 
law of thermodynamics, relating changes in the entropy to changes in the energy and other 
work terms (when appropriate). One attempts to understand this peculiar result by noting 
that horizons act as one way screens, causing information from a portion of the spacetime 
to become inaccessible to an observer on one side of it [Q. For example, an exterior observer 
cannot "see" what is going on within the event horizon of a black hole and a Rindler 
observer will never have access to the portion of Minkowski spacetime on the "other side" 
of her horizon. Notwithstanding this insight, the result itself remains highly non-trivial and 
intriguing because the appearance of a thermodynamic description suggests the existence of 
spacetime microstructure . In other words, it is possible that our description of spacetime 
is actually a description of the aggregate behavior of some as yet undiscovered microscopic 
degrees of freedom. 

This has led some investigators to suggest that gravity and perhaps even a continuum 
"spacetime" are emergent phenomena |^-|^. However, fundamental to the description of 
spacetime, at least in the current approaches, is the metric, which is determined by the 
equations of Einstein's relativity or of some generalization thereof, such as the Lanczos- 



Lovelock models in higher dimensions In an emergent picture of spacetime the metric 
must arise as an effective description, that is, out of the statistical mechanics of the more 
fundamental degrees of freedom of which it is made. How does this happen and how are 
its dynamics as encapsulated in, say, Einstein's equations to be reconciled with such a 
description? This motivates us to take a closer look at the structural relationships between 
our descriptions of dynamical systems and of thermodynamic systems. 

In this paper we ask the simpler question of whether there is anything truly "dynamical" 
about thermodynamics. That is, we ask if it is possible to describe thermodynamic processes 
as trajectories in the thermodynamic phase space in a way that is completely analogous to 
the description of a dynamical system. We argue that every thermodynamic system can 
be reformulated as a Hamiltonian system, analogous to the one describing a point particle 
moving in a curved space. The curved space is the configuration space of the thermodynamic 
system and its local geometry is captured by a matrix of functions, which incorporate all 
the physical attributes of the substance. If the matrix of functions is invertible, a quadratic 
Lagrangian description can be obtained. This leads to a natural interpretation of the matrix 
of functions as a metric on the configuration space: it determines a "distance" between 
equilibrium thermodynamic states. 

More precisely, the first law of Thermodynamics for quasi-static processes can always be 
written in the form 

df = p^dq^^ (1) 
where / represents a thermodynamic potential or the entropy and represent n intensive 
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or extensive variables of the system (including, possibly, the entropy), while each conjugate 
generalized force, 

P, = d,f, (2) 

is extensive if g'^ is intensive and vice-versa. As a simple example, consider the first law for 
an ideal gas in the form 

du = Tds — pdv (3) 

where u is the specific internal energy, s the specific entropy, v the specific volume of the 
gas and T,p are its temperature and pressure respectively. In this description the internal 
energy has the privileged role of / in (|l|) but an equivalent form 

du p ^ 

ds = — + ^dv (4) 

sees the specific entropy s playing that role. Therefore there is nothing special about / and 
the thermodynamic phase space is naturally odd dimensional. The first law can be viewed as 
the statement that physical thermodynamic processes may only occur on the hypersurface 
determined by the vanishing of the contact form 

u ■.= df - Pi^dq^. (5) 



Rajeev [|rT| has noted that there are many systems of variables in which the contact form will 
assume its canonical expression, as above, and all of them will be determined by requiring 
that 

aif,p,q)idf-p^dqn = df - P.dQ'^ (6) 

for any arbitrary function a on the 2n + l dimensional thermodynamic phase space. Trans- 
formations that take {f,p,q) to {f',P,Q) are the Legendre transformations of thermody- 
namics. Familiar and useful examples of Legendre transformations are transformations to 
the enthalpy, the Helmholz free energy and the Gibbs free energy. 

Here we will consider a restricted set of transformations, namely transformations of the 
2n dimensional subspace {q^,p^) — {Q^,P^) that leave the right hand side of (|l]) form 
invariant up to the addition of an exact form, i.e., 

df = P^dq'^ = Pf^dQ^ + da{p, q) (7) 

These are canonical transformations, under which f ^ f = f — cr. In section II we will 
determine a general equation for the generator of these transformations. Because they 
are canonical, the characteristic curves of the generating function will satisfy Hamilton's 
equations. In section HI, we confine ourselves to a special subset of these transformations 
for which an explicit form of the generating function can be obtained. This solution is given 
in terms of a matrix of functions. If the matrix is invertible it behaves like a metric on the 
configuration space and we show that there is a Lagrangian description for the evolution. 
This description is analogous to the Lagrangian description of a particle moving in a curved 
space defined by the metric. We argue that suitable metrics may be constructed starting 
from a complete set of equations of state, which in turn may be experimentally determined 
but are usually derived from a statistical model. Some examples are worked out in section 
IV and we end with a brief discussion of our aims and conclusions in section V. 
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II. CANONICAL TRANSFORMATIONS FOR THERMODYNAMICS 



Let U he a. local patch with coordinates {q^iP^) of a 2n dimensional manifold Ai on 
which is defined the contact structure (^. Consider a structure preserving transformation 
of the 2n coordinates, i.e., consider pa — ^ Pa{(l'^,Pfj.), (f — > Q^'isf'iPi^ such that 

df = Pt^dq^ = Pf^dQ^ + da{q, p) (8) 
where cr(g,p) is an arbitrary function of the coordinates, which gives 

0^P.§^ + |^ (9) 

op^ dpf. 

Transformations that satisfy these conditions are canonical transformations. An infinitesimal 
canonical transformation can be written as 

P,^P,=P^ + 5r7^^, qf^^Q^^ = q>^ + 6r6^ (10) 

where ri^{q,p) and e'^{q,p) are functions of and r is some parameter. Let (7{q,p) = 

6t \{q,p), then applying the conditions @ it is easy to see that 

de'' dX 



V/^ = -Pa 



dqf^ dq 



n /-.-.N 

Paj^ + = 11 

The solution of these equations can be determined in terms of the function F = Pa£°' + A as 

dF „ OF 

SO that F is the generating function of the infinitesimal canonical transformations. Finite 
transformations may be recovered by composing such infinitesimal transformations, that is 
by determining the integral curves of the vector field 

Thus the generating function F defines a one parameter family of curves. These are its 
characteristic curves and, in terms of the mock "time" parameter, r, introduced in (|10|), 
they satisfy the ordinary differential equations 

where the over dot represents a derivative with respect to r. We have thus recovered pre- 
cisely Hamilton's equations giving the classical trajectories of a particle described by the 
"Hamiltonian" F, 

dF 

F = P,g^ + A=p,— + A. (15) 
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If A = then / does not change during the process and F is extensive in the "momenta" , 

F(g^/ip^)=/iF(g^p^), (16) 

but if A 7^ 0, then 

f^f' = f-5TX (17) 
so that the transformations induce a change in / according to 

/ = -A. (18) 

Now the physical hypersurface of Ai is the one on which the contact form vanishes, i.e., 

dF 

f =P,r = P,Q^ = -^ + F (19) 



by (|15D and so it is determined by the condition that F = 0. 

III. A RESTRICTED CLASS OF TRANSFORMATIONS 

Let us henceforth consider transformations for which A = A(g). In that case, the following 
is a formal solution of (P!3|) for F, 



^(7^-(g)p^p, + A(g), (20) 



where g^'^{q) is an n x n dimensional matrix of arbitrary functions on the configuration 
space. The constraint, F = 0, defining the physical hypersurface translates into 



^jg^-{q)v,.Vv = -A(g) (21) 
and the evolution equations for this system are found from (|l^ to be 



BF 1 

P. = = -2^"'(^)^?"' - (22) 

where use has been made of (0). If in addition the matrix of functions g^^{q) is invertible, 
then it may be thought of as a metric on the configuration space of the thermodynamic 



system, for one finds that the constraint in (21) turns into the condition 



g^u^it = 1 (23) 

and this gives a natural distance 

dr^ = g^udq^dq" (24) 

between equilibrium states of the system. Furthermore, because r is an arbitrary parameter, 
the right hand side of the above equation is invariant under coordinate transformations 
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implying that the metric g^i, transforms as a covariant tensor of rank two. Using the first 
of (^), one has 

P,. = -mgA^W (25) 

and the second equation simphfies to 



X-'g^-'T^M. = -d,X (26) 



where "r^^^" is the Christoffel connection defining parallel transport induced by the metric 



'g. These equations can be recovered by extremizing the reparametrization invariant action 
S = J rfr£(g(r),g(r),r) = -J drXiq)^ g^,{q)qf^{T)r{r) (27) 
where the momentum is defined in the usual way as 

P. - §■ (28) 

The action (p7|) determines the minimum length path between the initial and final states in 
the thermodynamic configuration space with the conformal metric 

9f.u = X\q)g^,. (29) 

The mock time parameter, r, can be related to / by choosing A. Taking A = Aq (constant) 
forces r to be proportional to / according to (p^ . In this gauge, the action describing the 
evolution of our thermodynamic system becomes 

5 = -Ao j ^JgM dq^" dq^ = -Aq J dr, (30) 

which is precisely the action governing the evolution of a particle of mass Aq in a curved 
space defined by the metric 'g. 

The metric necessarily contains all the information about the substance, so it must be 
determined from appropriate additional considerations. Once determined, all processes are 
described by its geodesies. One possible choice of metric, that we will not adhere to in the 
examples of the next section, is obtained by noting that the quadratic form in (^) must be 
required to be positive definite. This is guaranteed by the second law of thermodynamics, 
i.e., by the concavity of the entropy function, s, if we take 



dq^^dq-- ^^^^ 



This metric has been proposed by Ruppeiner|T^ as it arises naturally out of classical ther- 
modynamic fluctuation theory. However, it requires a knowledge of the entropy as a function 
of the extensive variables, which is tantamount to a complete knowledge of the thermody- 
namics. Likewise, for the Weinhold metric which requires a knowledge of the internal 
energy. 
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From a thermodynamic point of view, a substance is characterized by a set of equations 
of state that are experimentally determined. For example, an ideal gas is completely char- 
acterized by the law of Boyle, Charles and Gay-Lussac together with a relationship between 
its internal energy and its temperature. These relations, although historically experimental, 
are also obtained directly from statistical models. In the statistical approach, a theoretical 
model is constructed and a partition function, Z{I3, g*), where (3 is the inverse temperature, 
is obtained, from which the internal energy and conjugate forces are determined according 

These n equations of state suffice to construct a family of invertible metrics all of which 
define a quadratic form of the kind g^'^PfiPu, whose value is a constant on the physical 
hypersurface. More precisely, suppose that we have n independent equations of state in the 
form 

fMP,=Pn (33) 

where pn are constant. The matrix / must be invertible so that it is possible to recover the 
momenta from the equations of state, 

P, = U-'rM)Pn. (34) 

If we take the p„ to define an orthogonal basis in a real vector bundle over the configuration 
space, with an invertible, constant matrix r/mn, then can be thought of as a vielbein, and 
one has the natural metric 

= r/-"/^/: (35) 
on M.. Furthermore, on the physical hypersurface, i.e., when F = 0, 



A = -^g^'^'P^Pu = -Vv"'''PmPn = Ao (36) 

is a constant determined by ff and p„ (we shall henceforth drop the subscript "0"). Taking 
rj to be diagonal, we find that g'^" is actually an n parameter family of metrics. As we will 
see in the following examples, they serve to parameterize the solutions. 

We assume that a complete characterization of the substance by means of n independent 
equations of state of the above form is available. Using these equations as our starting 
point, we assemble a metric according to (|35D, which automatically gives a constant A on 
the physical hypersurface. Below we illustrate the formalism for some common systems. 



IV. EXAMPLES 

Because the experimentally determined equations of state do not generally include the 
entropy, we single it out by letting / = s in the following examples. This is not necessary, 
but it is convenient for our purposes. Further, all of our examples are two dimensional and 
we will have two independent equations of state, so we take 

Vn^n = ("q ) (37) 

for arbitrary a and /3. 
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A. The ideal gas 



Take f = s, the specific entropy of tlie gas, and tlie first law in the form 



1 p 
ds = —du + —dv 
T T 



(38) 



In this representation, the coordinates are (u, w), the momenta Pu = '^/T, p^ = p/T and the 
ideal gas can be characterized by the two equations of state 



PuU 



PvV 



2 
k 



(39) 



where k is Boltzmann's constant and g is the number of degrees of freedom per molecule. 
According to (p^) this gives the two parameter family of metrics 



4 

n '"'^ 



(40) 



and, from the on-shell condition. 



A = 



I g'^k'^ k"^ 
4^ ^ ^ 



(41) 



We find the evolution equations 



u 



PuU" 



Pu 



pIu 



Pv 



PvV 

' I3^X'' 

pIv 

' /32A' 



(42) 



which will be seen to directly reproduce the equations of state and, furthermore, using the 
fact that s = — Ar, give the solutions 



u = Uq exp 



2a2A2 



V = vq exp 



ks 



(43) 



together with corresponding solutions for p^ and py, which follow from p9|). They can be 
inverted and the constants A, a and (3 eliminated using to give 



kin 



U \ V 



(44) 



The solutions in (|43|) describe a two parameter family of ideal gas processes. For instance 
the choice a'^ = gk/2, = k makes (1^) equal to the Ruppeiner (entropy) metric and 
describes a constant pressure process. A constant temperature process can be recovered in 
the limit as a — )■ oo and a constant volume process in the limit as /9 — t- oo. It is worth 
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noting that these trajectories are also recovered as geodesies of the metric in (|40|), or from 
the action 

^ = -A y dr^x^ + y^, (45) 

where x = alnu, y = (3\yiv and + = 1, as indicated earher, because the metric is flat. 
This is a consequence of the coordinate invariance of the formalism. 



B. The Van der Waals gas 



The treatment of the Van der Waals gas in similar, although the evolution equations 
appear more complicated at first sight. In the representation above, the Van der Waals gas 
is defined by the equations of state 



so fpSf) turns into 



\ 

and on the physical hypersurface, F = 0, 



a\ gk 
(^^ + ^) {v-h) = k 

a(v-bf {v-hf 



With this we find the following evolution equations 



u 



Pu 



Pv 



apl ( a\ 2apu ( apu 
M + - I - \ Pv^ — 



(46) 



(47) 



(48) 



Pv + ^] (v-b) 



(49) 



Now, employing the middle two equations, the first and last can be put in the form 



d / a 
dr \ V 



Pu ( , ay 
a^A V V/ 



d_ 

dr 



(pv + 



apu 



Pv + 



apu 



{v-b) 



(50) 
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from which it is clear that 



is recovered and, furthermore, the solutions 



u 



a a 
h Mo H I exp 



ks 



gks 



h+ {vq- h) exp 



(51) 



which may be inverted to obtain 

s = A; In 



^ 5/2 



^ DO 



(52) 



once use is made of (^81). As before, ( pTD can be recovered as geodesies of the metric (^7[), 
or from the action in (|45|) with 



X = aln {u ^ — j , y = (3 ln(t> — b). 



(53) 



C. Paramagnetism 

For definiteness consider spin ^ paramagnetism for which the first law will read 

ds = —dh — ^dm (54) 

where h = u + mh is the specific enthalpy, m is the magnetization and h is the magnitude of 
an applied, external magnetic field. In this representation the configuration space is made 
of the pair {h,m) and the conjugate forces are respectively ph = 1/T and pm = —b/T. A 
paramagnetic material does not have interactions between its dipoles and its enthalpy is 
completely independent of the magnetization, depending only on T. The simplest model 
of this non magnetic contribution would be to imagine that the paramagnetic molecules 
are oscillators, oscillating about their equilibrium positions. Therefore, if we suppose the 
enthalpy to be given by the law of equipartition, the material is completely characterized 
by the relations 



m 



1 A*^ 
/itanh — 



hph = gk 



(55) 



where fi is the magnetic moment of the particles. We can write these relations in terms of 
the phase space variables in the form 



k 1 -1 / ^\ , , 
Pm = tanh — , hph = gk 



They suggest the family of metrics 







_ I a2 

^ ' j^[Arctanh(m/;u)]-2 



(56) 



(57) 
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for which 



/ g2f^2 



X 



a" 



+ 



A;2 



on the physical hypersurface, F = 0. We find the equations 



h 



m 



Ph 



Pr. 



Pn 



/32A[Arctanh(m/ yu)]^ 

2 
h 



a^X 



Pm 



from which follow 



/3^Xfi{l — m2//i2)[Arctanh(m/yu)]^ 
and the solutions 



h = ho exp 



m 

— m 

2/i 



1 + m/ 



m/ /i 



+ - In 
2 



m 



/32;x2A2 



+ So. 



(58) 



(59) 



(60) 



The second determines the magnetic contribution to the specific entropy. That the entropy 
is the sum of the two contributions follows by simply eliminating the arbitrary constants, A, 
a and (3 in the two equations above, using (^8|). Once again, these solutions are geodesies 
of (0). They can also be obtained from the action in ( ^5]) with 



x = a\nh, ?/ = /3|mtanh ^ (m//i) + ^ln[l — (m//i)^]| 



(61) 



D. Kerr Black Hole 



The Kerr black solution is one of four black hole solutions in general relativity and 
describes a neutral, rotating black hole. It is completely characterized by its axisymmetry 
and by two parameters, viz., the mass of the black hole, M, and its angular momentum, J. 
It is described by the spacetime metric (we take c = 1 and G = 1) 



r^r 



l-^]dt 



2r,ra sin^ 9 



dtd^-^dr^-p^de^- (r^ + a^ + 



sin^ 9 difi"^ (62) 



P / P 
where is the Schwarzschild radius, = 2M, and 

a = J/M 

p"^ = r^ + cos^ 9 

A = — + 

There are two surfaces of interest: an inner surface (the event horizon) occurring at 

rh = ^{rs + V'r2 - 4a2) 



(63) 



(64) 



11 



and an outer surface of infinite redshift at 

Te = 2{'''^s + v^r^ — Aa? cos^ d) (65) 

The two horizons meet at ^ = and the region between them is called the ergosphere. 
Within the ergosphere, a test particle must co-rotate with the mass M , with the angular 
velocity 

n = (66) 

A quantum field placed in this background is well known to acquire a temperature. This 
temperature is proportional to the acceleration of the null Killing vector on the horizon, 
known as the surface gravity, k, of the hole, 

T= — = (67) 

The first law of black hole thermodynamics can be written in the form 

dS = ^dM - ^dJ. (68) 

From the thermodynamic point of view, the configuration space is two dimensional, spanned 
by the Arnowitt-Deser-Misner mass, which plays the role of the internal energy, and the 
angular momentum. The corresponding conjugate forces are pm = 1/T and pj = —Q/T. 
Although these are natural variables for the system, it is convenient to transform to the 
configuration space (r/j,a), writing the first law of black hole thermodynamics in terms of 
these variables as 

ds = Phdrh + Pada (69) 
instead. The equations of state, (|66|) and (j§7h, are then equivalent to the statements that 







Ph = 2nrh and Pa = 2na and one finds 



9'" = I °? 1 I , (70) 

which yields 



a2 /32 

on the physical hypersurface. The evolution equations 

Ph 



^ = -2^A/i + 4- (71) 



Pa 



Ph 



pI 
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(72) 



have the solutions 



Ph Pa 2 47rS 2 2 ^TTS 

- = const., - = const., = + a = ^ 



and the black hole entropy 

s = n{rl + a^) + sq (74) 

is recovered as the simple sum of contributions from and by eliminating the constants 
using (^). As in our previous examples, these solutions are geodesies of ( [TOl) and can also 
be obtained from the action in (|^ with 

x = a^rl y = f3^a^ (75) 



V. DISCUSSION 

If spacetime is indeed emergent, its microscopic degrees of freedom quite possibly live 
at scales on the order of the Planck length. Therefore, not only do we presently have no 
experimental access to to them but we are most likely never to have it. It would seem 
that the best we could hope to empirically justify is an ever more precise description of 
the thermodynamics of this microstructure and, in turn, such a sharpened description may 
eventually lead to a better understanding of the fundamental constituents of spacetime. 
Progress along these lines can be made only once we have a clearer picture of the connections 
between our current geometric description of spacetime and thermodynamics. In this paper 
we have begun to address this issue by taking a closer look at thermodynamics. 

We were able to show that at least a subset of the possible transformations on the 
thermodynamic phase space lead to the description of thermodynamics processes as geodesies 
of a family of metrics defined by the equations of state of the substance. In the entropy 
representation, it is the entropy that serves as thermodynamic "time" . In general, it is the 
"preferred" function / that plays this role. Although all the metrics in the examples we 
have considered are flat there is no reason to expect that this is universally so, particularly 
in higher dimensional systems, eg., systems with variable contents or charged and rotating 
black holes. It would be of considerable interest to study the non-trivial geometry of these 
systems and correlate their geometric properties with their thermodynamic behavior. Such 
work, based on the Ruppeiner or Weinhold metric, has been attempted [|l4l-[l6| and our work 
can be seen as providing additional motivation for it. 

It is interesting that the family of relevant metrics is completely recovered simply from the 
equations of state. This is in fact what one should expect in an emergent picture of spacetime: 
that the geometric formulation is simply a way of specifying the spacetime equations of 
state. We hope that, with more work, this insight may help to "design" geometric models of 
emergent gravity that are better behaved in the ultraviolet and perhaps even to eventually 
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understand spacetime's microstructure better in an information theoretic way. 
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